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Precision computation of hadronic physics with lattice QCD is becoming feasible. The last decade 
has seen percent-level calculations of many simple properties of mesons, and the last few years 
'yr' have seen calculations of baryon masses, including the nucleon mass, accurate to a few percent. 

As computational power increases and algorithms advance, the precise calculation of a variety of 
' more demanding hadronic properties will become realistic. With this in mind, I discuss the current 

JlJ lattice QCD calculations of generalized parton distributions with an emphasis on the prospects 

for well-controlled calculations for these observables as well. I will do this by way of several 
examples: the pion and nucleon form factors and moments of the nucleon parton and generalized- 
T— I parton distributions. 
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1. Introduction 

Over the last decade, lattice QCD has proven itself capable of percent-level calculations. This 
is a milestone that establishes lattice field theory as a powerful tool for performing reliable and 
quantitative computations of nonperturbative QCD phenomena. However, it is important to un- 
derstand that the computational demands required for such calculations depend strongly on the 
quantities of interest. Consequently, the most impressive lattice calculations to date have been 
Umited to arguably the simplest QCD observables, mostly meson properties. Recently, there has 
been progress on the determination of the baryon spectrum. In particular, the nucleon mass can 
now be calculated nonperturbatively from QCD at the few-percent level, raising the prospects for 
calculations of more challenging observables that have long been sought from lattice QCD. 

Even with the successful computation of the nucleon mass, the calculation of nucleon form 
factors remains a challenge for lattice QCD. Significant progress has been made in the last few 
years. In particular, recent computations with nearly physical pion masses represent a major break- 
through for lattice QCD and have been essential to the successful determination of the nucleon 
mass, which is clearly a necessary first step. However, there are additional sources of uncertainty 
that occur for calculations of nucleon matrix elements. With nearly physical pion masses becoming 
more common, the focus is turning toward fully controlling all uncertainties relevant to the compu- 
tation of matrix elements in order to perform calculations of form factors that can be quantitatively 
compared to experimental measurements. 

In these proceedings, I start by briefly reviewing several examples of high-precision lattice 
calculations and the recent precision determinations of the masses of baryons, in particular the 
nucleon. I then turn to the pion form factor as a simple example of the more challenging quantities 
that our community is interested in. The calculation of this observable provides a glimpse of the 
progress that we can hopefully expect for nucleon form factors in the years to come. However, 
nucleon calculations are still an open issue, so I finish by looking at just a few representative 
examples of the generalized parton distributions of the nucleon. 

2. High-precision lattice QCD calculations 

There is now a small but growing list of quantities that can be calculated using lattice QCD 
with percent-level accuracy. Figure 1 shows several. In order to understand if and when such 
precision will be brought to bear on more challenging observables, it is important to note that the 
earliest calculations in figure 1 were only performed within the last decade. Percent-level precision 
has been reached for additional quantities since then, but almost all of these are related to simple 
properties of stable mesons. Despite these limitations, the results in figure 1 represent a signif- 
icant accomplishment of the lattice community and demonstrate that high-precision calculations 
ai^e possible using lattice QCD. 

3. Precision lattice QCD computation of baryon masses 

There is a notable absence of baryon properties in figure 1. For a variety of reasons, calcula- 
tions involving baryons require significantly more computational resources that those for mesons. 
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Figure 1: Examples of high-precision lattice QCD 
calculations. Several observables have been calcu- 
lated with a total precision of one percent or bet- 
ter, illustrating that lattice QCD is capable of high- 
precision computations. The quantities shown here 
are averages of several lattice calculations [1,2]. 




Figure 2: An example precision computation of 
baryon masses. Lattice computations involving 
baryons are more challenging than mesons and have 
only recently reached the few-percent level. The 
masses from BMW [3] are shown, but PACS-CS [4] 
and ETMC [5] find similar results. 



In particular, nucleon properties appear to have a more significant dependence on the pion masses 
used for their determinations. This represents a computational threshold that afflicts the calculation 
of nucleon observables more than others. Nonetheless, calculations completed in just the last few 
years were the first to demonstrate some quantitative control over the ground-state baryon masses 
shown in figure 2. The determination of the nucleon mass niN at the few-percent level is an impor- 
tant advance forward for the lattice community. It raises the prospects that the more challenging 
task of computing nucleon matrix elements will be within the reach of lattice QCD. However, we 
should bear in mind that the calculation of m/^ has only been possible in the last few years. This 
suggests that nucleon matrix elements may very well require yet more time. 



4. Lattice QCD calculation of the pion form factor 

The pion electromagnetic form factor has recently started to yield to quantitatively controlled 
lattice calculations. It serves as an example of what to expect for nucleon form factors in the years 
to come. Figure 3 shows a summary plot of lattice QCD calculations of the charge radius of the 
pion {r^)n, defined by 



{K,p\J^\K,k) = {p^+k^)F{Q^) Q^ = -{p-kf l/)n = -^ 



e-=o 



where = Y^fQfjli, jji = qfY^Qf and Qf is the electric charge for each quark flavor /. There ap- 
pears to be no agreement amongst the lattice calculations, some of which do and some do not agree 
with the experimentally measured value of {r^)^. As I will argue shortly, this plot is misleading, 
but I have chosen to examine it as an example of how lattice QCD calculations progress. 

To understand what is happening in figure 3, we can crudely plot the results for {r^)n versus 
their publication date. This is done in figure 4, which shows a time dependence in the calculations 
with later results agreeing with the measured {r^)^- More importantly, after selecting those cal- 
culations that have accounted in some way for the relevant sources of uncertainty (highlighted in 
figure 4), we see a consistent picture emerge. All well-controlled lattice calculations agree with 
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Figure 3: An example of a typical but mislead- 
ing summary plot. This summary of the pion 
charge radius {r^)K, similar to the one in [6], il- 
lustrates the dangers in comparing all lattice com- 
putations without regard to controlled uncertain- 
ties. The values are from PDG [7], ETMC [6], 
JLQCD [8], JLQCD/TWQCD [9], LHPC [10], 
QCDSF/UKQCD [11] and RBC/UKQCD [12]. 
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Figure 4: Time history of lattice calculations of 
the pion charge radius {r'^)K- This time history il- 
lustrates that lattice QCD calculations with well- 
controlled uncertainties represent definitive predic- 
tions of QCD that can be safely compared to exper- 
imental measurements. The impression given by 
this plot should be contrasted with that given by 
figure 3, which fails to distinguish those calcula- 
tions with fully controlled uncertainties. The val- 
ues are the same as in figure 3 and the two high- 
lighted calculations are, starting with the earliest, 
by QCDSF/UKQCD [11] and ETMC [6]. 



each other, as should always happen, and, in this case, agree with the experimental measurement. 
At this point we are left with only two calculations that meet this burden [6, 11], but several of the 
calculations in figure 4 are ongoing and there are new calculations [13] underway. We should ex- 
pect in the next few years to have multiple precise and reliable lattice calculations of {r^)^- As with 
the results in figures 1 and 2, the easiest observables are calculated first, but this is to be expected 
from large-scale numerical calculations that face multiple computational thresholds that depend on 
the quantities of interest. The calculation of (r^);r is not only interesting in its own right but also is a 
good example of how nucleon form factor calculations may progress as the computational barriers 
are crossed for those observables as well. 

IVIoving beyond the extreme 2^ — limit characterized by the slope of the form factor, lattice 
calculations are also exploring the non-zero but still low regime of < 1 GeV^. Various 
uncertainties become more challenging with increasing Q^, but the low behavior of F{Q^) can 
also be calculated with relatively well-controlled uncertainties. This is shown in figure 5 along with 
the experimental results in the low region. The agreement between lattice QCD and experiment 
is again rather compelling and bodes well for the eventual determination of nucleon form factors. 

5. Status of nucleon form factors from lattice QCD 

Lattice calculations with nucleons remain a challenge due to the extra computational demands 
that are required. As such, there is not yet a compelling success story for lattice determinations of 
nucleon structure, but this needs to be understood in context. It has only been in the last few years 
that the lattice community has developed sufficient algorithms and garnered enough computing 



4 



Form factors from lattice QCD 



Dm B. Renner 



0.8 



0.6 



a 



0.4 



0.2 




o CERN 
□ DESY 
A JLAB 

■ lattice QCD (ETMC) 



0.2 



0.4 0.6 
Q' [GeVI 



0.6 
0.5 



A 0.4 



A 0.3 
0.2 
0.1 




• CODATA 
o ^i-H 

— HBxPT LO 
■ lattice QCD (QCDSF) 

* lattice QCD (LHPC) 



finite-size effects for small m 



■■5 



0.05 



0.1 0.15 
[GeV^] 



0.2 



0.25 



Figure 5: Comparison of the pion form fac- 
tor F{Q^) from experiment and lattice QCD. 
ETMC [6] has performed a lattice calculation with 
reasonably controlled uncertainties of the low 
dependence of F{Q^). This is an example of the 
level of control sought for the more challenging 
nucleon form factors. The experimental measure- 
ments are from CERN [14], DESY [15, 16] (reana- 
lyzed in [17, 18]) and JLAB [17]. 



Figure 6: The proton isovector charge radius 
The measured values use CODATA [19] or 



li-H [20] for {r^)p with PDG [7] for (r^)^, to form 
('"^)p = {>'^)p^ {''^)n- These are then matched to 
chiral perturbation theory [21]. The lattice results 
are from QCDSF [22] and LHPC [23]. Requiring 
well-controlled uncertainties likely eliminates the 
apparent tension between lattice QCD, chiral per- 
turbation theory [21] and the measurements. 



power to calculate even m/v at the few-percent level. Thus it should not be surprising that current 
calculations of nucleon structure have produced mixed results, with some encouraging agreements 
and some noticeable disagreements. As illustrated by the pion form factor, this is how we expect 
lattice calculations to look before fully controlled uncertainties are achieved. 

With this in mind, I simply summarize the current status of nucleon calculations with a few 
examples. To compare with the discussion of the pion charge radius, I start with the proton charge 
radius defined similarly. 



{P,p\4\P,k) = KlF({Q') +KlFl{Q') {r% = ^""'^^'^ 



2^=0 



where K^^ are kinematic functions of p and k. From the outset, I have to make a sacrifice and instead 
focus on the isovector radius (r^)^^'^ = — {^'^Yp- This restriction eliminates so-called discon- 
nected diagrams that are an additional source of uncertainty. Furthermore, I avoid any statements 
about extrapolations to the physical point and simply examine the lattice results as functions of the 
artificially heavy pion mass m-j^ used in the calculations. These are shown in figure 6. To avoid 
unnecessary clutter, only two lattice calculations are shown as examples, and only the uncertainties 
due to numerically integrating the path integral are shown. Considering all results in figure 6 with- 
out regard to the remaining uncertainties in the problem, there would be an apparent disagreement 
between lattice calculations and the measured value. Additionally, there would seem to be a signif- 
icant disagreement with the expectations of heavy baryon chiral perturbation theory. However, as 
with the early calculations for the pion form factor, the failure to accurately control all uncertainties 
in the lattice calculation may very well be the explanation for both of these observations. 

A detailed discussion of the uncertainties for (r^)p^^ is beyond the scope of these proceedings, 
but I will examine one uncertainty that may be a concern: finite-size effects. Lattice calculations 
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Figure 7: The isovector momentum fraction 
{xYp-''. The lattice results are QCDSF [24] and 
ETMC [25]. The results from global analyses are 
from [26] using [27-32]. The leading -order results 
from chiral perturbation theory [33,34] are matched 
to the largest and smallest results from the global 
analyses to indicate the impact of the variations of 
these results. The red arrows indicate the possible 
impact of the uncertainty, suggested by ETMC [35] 
and LHPC [36], due to excited-state pollution in the 
lattice calculations as described in the text. 




Q [GeV-] 

Figure 8: The isovector A'^^^''{Q^) generalized 
form factor. The moments of the generalized- 
parton distributions of the proton can be calcu- 
lated using lattice QCD techniques. The general- 
ized form factor A2q{Q^) shown here is related to 
the first x-moment of the sum of the //(x, ^ , A^) and 
E{x,^,A^) GPDs in which = A^. The results 
from [37] are shown as an example. The absence of 
experimental results emphasizes the prospects for 
genuine predictions from lattice QCD once all un- 
certainties are reliably controlled. 



are performed at a finite piiysical volume Lr" and correspondingly the allowed momentum modes 
are discretized. Forming the derivative of the form factors requires taking finite differences in Q^. 
Furthermore, this quantity is expected to diverge in the chiral limit, so we may reasonably expect an 
enhanced sensitivity to the volume, which acts as the infrared regulator for lattice QCD. Due to the 
demanding nature of nucleon calculations, the large volume limit is difficult to study, especially as 
nin is decreased. Consequently, there is not much evidence for or against a strong finite-size effect 
for (r^)^p'^. In lieu of a full study of the L dependence, we can require the usual rule of thumb 
mjtL > 4. This restriction in figure 6 clearly alleviates the disagreement with the measured value 
for {r^Yp^'^ and substantially weakens the tension with the expectations from chiral perturbation 
theory, suggesting that calculations with larger L may help resolve these puzzles. 

Nucleon form factors are challenging enough; turning to moments of parton and generalized- 
parton distributions increases the difficulties further. In figure 7, 1 show two example calculations 
of the average longitudinal-momentum fraction x in the proton (x)^^^ given by 

where the K'^y are kinematic functions, O^p^y = qfY{pDyjqf are symmetric traceless twist-two op- 
erators, and q^{x) can be related to the quark (anti-quark) parton distribution functions (PDFs) for 
X > (x < 0). Generalizations to higher moments {x")p are possible. The isovector combination is 
used for the same reasons as discussed regarding {r^Yp '^- There is an apparent disagreement with 
the measured value, but again the issue boils down to fully controlling the relevant uncertainties. 
An important uncertainty for (x)^^'^ appears to be a finite-size effect due to excited-state pollution. 
(This is a technical issue related to having a large enough separation between correlation functions 
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in the calculation.) In the case of (x)^-^ there are indications from lattice calculations that this 
effect is occurring. The red arrows in figure 7 show the sort of corrections suggested by direct 
calculations. These corrections are not universal and can not readily be applied to the results of 
other computations, so we can only note that the sign of the correction and its the m^i dependence 
seems to reduce the disagreement with the measured (x)^^'^'. 

I have focused on {r^)p^'^ and (x)^^'^' because they are commonly used as benchmarks for 
lattice calculations of nucleon structure. Just as the calculation of mm was an essential achievement 
for making progress, the successful determination of {r^Yp'^'^ ^^'^ other similar observables 

is a necessary step towards the broader structure program envisioned by the lattice community. As 
just one example of that effort, I consider the generalized-parton distributions (GPDs). As for the 
PDFs, lattice calculations focus on moments in x. The first moments of the H and E GPDs are 
given by the generalized form factors A2o( 2^), ^20 (2^) and €2(0^) ■ One example is 

where the definitions of the H and E GPDs and further details are available in [38] . There have been 
several lattice calculations of these form factors, and one recent example is shown in figure 8. The 
calculation of '^(2^) is an extension of (x)^^^ and all issues relevant for (x)^^'^ are expected to 
occur for ^20^(2^) as well. In fact, A2o'^'(0) = (x)^^''. One goal of the long-term lattice effort on 
nucleon structure is the determination of the low moments of all the nucleon GPDs. The restriction 
to moments in x is a significant limitation to PDFs, but for GPDs the additional information on 
the = A^ dependence adds valuable information on the GPDs that is complimentary to that 
which is accessible from experimental measurements. Thus well-controlled calculations of the 
dependence of form factors like A2o(2^) will yield genuine predictions from lattice QCD for the 
nucleon and will open several avenues to rich physics topics including the spin decomposition and 
ti^ansverse structure of the nucleon. 

6. Conclusions 

In the last decade, lattice QCD has shown itself to be capable of precision calculations. The 
initial successes were understandably limited to those observables for which all uncertainties could 
be controlled with the least computational resources. However, recent calculations, such as the 
determination of m^, indicate that lattice computations of more demanding hadronic quantities 
should become feasible in the years to come. 

The pion form factor is currently a well-determined quantity from lattice QCD, with further 
improvements expected. It stands as an example of the way forward for equally well-controlled lat- 
tice calculations of nucleon structure. There will be continued emphasis on benchmark observables, 
such as (r^)J^^'^' and (x)J^^'^', in order to establish control of all relevant uncertainties for calculations 
of nucleon matrix elements. This renewed focus on carefully accessing each uncertainty is already 
shedding some light on cun^ent puzzles facing lattice calculations of nucleon structure and will 
eventually lead to precise and reliable calculations that can safely be compared to experimental 
measurements. However, the real promise of lattice calculations of nucleon matrix elements ul- 
timately lies in determining observables that lie beyond the reach of measurements yet offer the 
possibility of insight into the deeper mechanisms behind nucleon structure. 
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